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Abstract 

In this paper we revisit the problem of a (non self- avoiding) polymer chain in a 
random medium which was previously investigated by Edwards and Muthuku- 
mar (EM) As noticed by Gates and Ball (CB) Q there is a discrepancy 
between the predictions of the replica calculation of EM and the expectation 
that in an infinite medium the quenched and annealed results should coincide 
(for a chain that is free to move) and a long polymer should always collapse. 
CB argued that only in a finite volume one might see a "localization tran- 
sition" (or crossover) from a stretched to a collapsed chain in three spatial 
dimensions. Here we carry out the replica calculation in the presence of an ad- 
ditional confining harmonic potential that mimics the effect of a finite volume. 
Using a variational scheme with five variational parameters we derive analyt- 
ically for d < 4 the result R {g |ln/i|)"^/(^~'^) ~ {g\nVy^l'^^~'^\ where R is 
the radius of gyration, g is the strength of the disorder, ^ is the spring con- 
stant associated with the confining potential and V is the associated effective 
volume of the system. Thus the EM result is recovered with their constant 
replaced by \nV as argued by CB. We see that in the strict infinite volume 
limit the polymer always collapses, but for finite volume a transition from a 
stretched to a collapsed form might be observed as a function of the strength 
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of the disorder. For d < 2 and for large F > F' ~ exp(c/2/(2-<i)i^(4-<i)/(2-d)) 
the annealed results are recovered and R ^ {Lgy/^'^^^\ where L is the length 
of the polymer. Hence the polymer also collapses in the large L limit. The 
1-step replica symmetry breaking solution is crucial for obtaining the above 
results. 

PACS number(s): 05.40-a, 75.10.Nr, 36.20.Ey, 64.60.Cn 
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I. INTRODUCTION 



There has been much interest in recent years in the properties of polymer chains in a 
quenched random environment ■ This problem is directly related to that of a quantum 
particle in a random medium and to that of a flux line in a type II superconductor in 
the presence of random columnar defects |^,^ as will be made clear below. Thus its general 
application makes it important for a variety of physical situations. 

The quantities of interest for the polymer problem are the free energy and the radius 
of gyration of a chain in a quenched white-noise potential. Here we consider only the case 
of a non self-avoiding chain. Gates and Ball |^ gave a beautiful intuitive argument to the 
effect that a Gaussian chain situated in an infinite random medium is always collapsed in the 
long-chain limit. Their argument goes as follows: Gonsider a white-noise random potential 
v{x) of zero mean whose probability distribution at each site is: 

P(f(x)) oc g'^/^ exp{-vy2g). (1.1) 

If we now coarse-grain the medium and denote by v the average value of the potential 
over some region of volume a, Then the coarse-grained potential will have the distribution 

Pa(v) OC (g/a)-'^' expi-avy2g). (1.2) 

Gonsider a polymer chain situated in the random potential, and assume that it shrinks 
into a volume a corresponding to a place where the mean potential v takes on a lower 
value than usual. In this situation the free energy of the chain is crudely estimated to be 
(neglecting all numerical factors): 

F{a,v) = L/R^ + Lv + av^/2g. (1.3) 

Here L is the length of the chain (number of monomers), R is the radius of gyration (or end 
to end distance) and the volume a is related to R via a = i?'^ in rf-spatial dimensions. The 
first term on the r.h.s. is an estimate of the free energy of a long chain confined to a region 
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of size R in the absence of an external potential (see e.g. |T0|, Eq. 1.12). The second term is 



just the potential energy of the chain in the random potential of strength v. The third term 
arises from the chance of incurring a random potential of strength v. The quantity InP(w) 
gives an associated effective entropy for the system. Minimizing this free energy over both v 
and a determines the lowest free energy configuration. Minimizing with respect to v yields 
V = —Lg/a. Substituting in F gives: 

This shows that for any d > 2, F — >■ —oo as R —>■ 0. Thus the mean size of the chain is zero, 
or in the presence of a cutoff of a size of one monomer, 

i? ~ 1, d>2. (1.5) 

For d < 2, the free energy has a minimum for 

(L^?) 1/(^-2) d<2, (1.6) 

which in the long chain limit (L > 1/g) cuts off again at i? ~ 1. These results are the 
same as those for the case of an annealed potential that is able to adjust locally to lower the 
free energy of the system. The reason is that for an infinite system containing a finite (even 
though long) chain, space can be divided into regions containing different realizations of the 
potential, and the chain can sample all of these to find an environment arbitrarily similar 
to that which would occur in the annealed situation. 

These results stand in contrast to the replica calculation of Edwards and Muthukumar 
(EM) [jl|, who found that for a long chain 

when (/^/(^"'^^L — > oo , whereas R ~ L^/^ when g'^/'^^~'^^L — > 0. Note that the result ( p77| ) is 
independent of L as opposed to Eq. (|T]^). To reconcile the two apparently different results. 
Gates and Ball argue that the quenched case is different from the annealed case only for 



the case when the medium has a finite volume V. In a finite box, arbitrarily deep potential 
minima are not present. Instead the most negative v averaged over a region of volume 
a <^V occupied by the chain, is approximately (keeping only leading terms in the volume 
V) given by solving the equation (the l.h.s. of which represents the area under the tail of 
the distribution) 

r dy Pa{y) c:^ ^ , (1.8) 

J —oo V 

which yields 



(1.9) 



This expression when plugged into Eq. (|L3|) leads to (Note that the last term in ( |1.3| ) just 
becomes a constant independent of R) 



When this free energy is minimized with respect to R it gives rise to 

i? ~ (^lnV)-i/(^-^\ d<4 (1.11) 

which agrees with Eq.( p..7D and also with simulations performed on a chain in a random 
medium of a fixed finite volume @]. However it is not clear from this explanation why the 
replica calculation which has been done for an infinite system |]l[ gives rise to the finite 
volume result. To shed light on this question we will show in this paper that the reason for 
the discrepancy is the fact that EM used a variational calculation which relies on a single 
variational parameter. We show specifically that the single parameter variational solution 
is inconsistent. 

What we will do first is, instead of considering a system in a finite volume which is 
hard to solve, introduce an external harmonic potential (with a spring constant fi). Such an 
attractive potential has the effect of confining the chain to a finite distance from the origin 
since the energy cost to wonder far away from the origin of the potential is high. A system 



in a harmonic potential is easier to solve than a system in a finite box. It also corresponds 
directly to the problem of a flux line in a type II superconductor where the cage potential 
felt by a flux line due to its neighboring flux lines can be modeled by a harmonic potential 
(see below). In addition, we introduce more variational parameters, three for the case of 
replica symmetric parametrization and flve for the case of replica symmetry breaking (RSB). 
These extra parameters have physical signiflcance as will be discussed below. We will then 
use the replica method and the variational approximation to tackle the problem and obtain 
the free energy and the radius of gyration. For flnite /i we flnd that R is independent of L 
(the chain length) and as the disordered strength is increased from zero, R is decreased from 
its initial //-dominated value according to the relation R ~ {g\\n ^\)~^^'^^~'^'^ (which agrees 



with Eq. ( |1.11|) since the effective volume available to a system in a harmonic potential is 
In y ~ |ln /i|). 

On the other hand, if we try to take the ultimate /i — ^ limit (which is the case 
originally studied by EM), the previous solution becomes invalid and the chain collapses 
for d > 2. For d < 2 the annealed results are obtained in the /i — > limit as given 



above in equations ( |1.4|jl.5| , prBD . This occurs speciflcally because of the extra variational 



parameters used beyond the single variational parameter used by EM. We also demonstrate 
the importance of RSB for obtaining the correct physical results. (The relevance of RSB to 
this model was recognized by Haronska and Vilgis 0, but unfortunately their calculation 
still predicted a constant coefficient of proportionality in the relation R ~ {cgY^^^^~'^\ that 
although differs from the EM result does not contain the correct InV^ dependence.) 

To define the model of a polymer chain in a random potential plus a fixed harmonic 
potential we use the Gaussian chain approximation to write: 



H 



/ du 




(dR{u)\ 


Jo 




[ du J 



[1.12) 



were R(u) is the ci-dimensional position vector of the chain at arc- length u {0 < u < L), n 
governs the strength of the harmonic potential and V^(R) is the random potential satisfying: 

{V{R)) = 0, {V{R)V{R'))=gS^'\R-R'). (1.13) 
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We can actually consider a wider class of random potential correlations characterized by a 
function /: 

{V{K)Vm) =gdf (^^^^^) , (1-14) 

where /() is some given function. In Eq. ( |1.12| ) we choose the units such that u is dimen- 
sionless and so L is the length of the polymer in units of the Khun bond step b. The "mass" 
M is inversely proportional to where f3 = l/ksT (in (i-dimensions l3M = d/h'^). The 
case R(0) = R(i^) corresponds to a closed chain. 

The partition sum is given by the functional integral 



Z(R,R',L,/5) = / I?R(m) exp(-/5ff). (1.15) 

JR(0)=R 

We further define a boundary-free partition sum (for a closed chain) by 

Z(L, 13)= JdR Z(R, R, L, (3), (1.16) 
and the free energy is given by 

f3F = -\nZ{L,(3). (1.17) 
The correlation function of interest is 

C{i) = l{{{R{i)-R{0)r))^ , (1.18) 

where 

1<£<L. (1.19) 



The first average in Eq. ( |1.18| ) is the thermal one with a Boltzmann weight exp(-/5if) and 



the second average is over the realizations of the random potential. For the range of i given 
by Eq. ( |1.19| ), the boundary conditions on the chain, e.g. open or closed are not important 
for the behavior of C{i). 

For the case of no disorder (i.e. g = 0) the correlation function is given by 



Co{i) = (l - exp(-£V^)) . (1.20) 

We see that in the hmit /i ^ 0, Co(^) ~ ^ / PM ~ (b'^/d) i, which corresponds to pure 
diffusion of the chain (random walk). From the relation 

we see that the polymer chain is confined to a volume of size V satisfying 

\nV^^\\nfi\, (1.22) 

for small fi. 

The mapping of this problem to a vortex line in an harmonic cage potential and random 
columnar defects is such that the arc-length u corresponds to the distance z along the c-axis 
(assuming this is also the direction of the magnetic field), M — e; = eo/7^, which is the line 
tension of the flux line and 7^ = mz/m± is the mass anisotropy. R is a two dimensional 
vector in the a — b plane of the superconductor The harmonic potential plays an 

essential role as a reasonable approximation to the cage potential that a vortex line feels 
due to the repulsion by its neighbors. Thus fi ~ eoi?/$o where B is the magnetic field and 
$0 is the fiuxoid. 

There is also a mapping into the problem of a quantum particle in a random potential 
+ a harmonic potential. This mapping reads 



P 1/h, L Ph, (1.23) 

and p(R, R',/5) = Z(R, R',L = ph, P = 1/h) becomes the density matrix of a quantum 
particle at inverse temperature p. The variable u represents the Trotter (imaginary) time. 
In this case M corresponds to the mass of the particle. 



II. THE VARIATIONAL CALCULATION 



In order to average over the quenched random potential we use the replica method. After 
introducing n-copies of the chain and averaging over the random potential one obtains 
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(2.1) 



with 



Hn =Io'du ELi[f (^) 

- f lo du /o^ du' Eab S^'^ {Ka{u)-Mu')) 



(2.2) 



Here we used the delta function potential (to make contact with EM) , but later we will show 
how to generalize to a general correlation. It is useful to replace the delta function by the 
equivalent expression: 

dk 



{Ka{u)-R,{u')) = J exp (^k ■ (R,(w)-R,(w'))) • 



For a general correlation (see Eq.( |1.14[ )) we can write 



(2.3) 



f[iRaiu)-Rbiu)y/d 
dk 



J dyf{y'/d) J ^exp(-zk-y) exp (zk- (R,(«)-R,(«'))) • 



(2.4) 



In order to proceed we use a quadratic variational Hamiltonian to be the best approxi- 
mation to Hn. This is given by 



lo du E:=i 



i Jq du Jq du' J2ab <lab{u - M')Ra(M) ■ Rfe(M'), 



(2.5) 



where qab{u) are n x n variational functions to be determined, with n — > at the end. The 
best variational Hamiltonian is determined by the stationarity of the variational free energy 
which is given by ||TT|Jl|Jl^: 



n{F)j^ = {Hn-hn)h^-^ln J PRi • ■ ■ PR, exp(-/3/i„). 



(2.6) 



The general equations satisfied by qab{u) where discussed in Refs. 1^,0,3. We showed that 
although the diagonal elements qaaiu) must depend on the arc-length variable m, the off 
diagonal elements Qa^b which are spin-glass like order parameters can to be chosen to be u- 



independent; in other words there is a consistent solution of the variational equations with 
these properties. (The existence of a time-persistent part to the off diagonal elements of 
Qab is well known in the investigation of quantum spin glass systems [|I^ and is crucial for 
the capture of the correct physics in such systems.) In Refs. [|],^ we proceeded to solve 
the equations approximately for the case of a non-zero confining harmonic potential char- 
acterized by a spring constant /i 7^ 0. For a quantum particle at not too low a temperature 
(equivalent for moderate values of L in the polymer problem) we obtained a numerical so- 
lution of the equations |^ for different types of correlations of the random potential. In 
Ref. we considered the limit of large L and finite /i (in the context of the vortex line 
problem), and for c? = 2, under certain approximations obtained an analytical solution to 
first order in g (the strength of the disorder). Here we would like to consider the whole range 
of disorder for large L and also investigate the limit /i — 0. Our goal is also to make contact 
with the calculation of EM. Hence we will start with a somewhat simpler approach with 
a finite number of variational parameters in lieu of the infinite number of such parameters 
introduced in our previous work. As will turn out this is appropriate for the current problem 
and allows us to solve everything analytically without any further approximations. 
EM considered only the case of /i = and chose 

qab{u-u') = —5ab5{u-u) (2.7) 

y 

where g is a single variational parameter. (In an Appendix they considered a slightly more 
general form but it is still proportional to 5{u — u') ). Here we claim that we need to 
introduce static (w-independent) off-diagonal elements for qat and also add a static diagonal 
part. This will help capture the correct physics of the problem as in the case of the quantum 
spin glass systems mentioned above. Thus we chose: 

qabiu - u') = -6ab ((A - /i) 6{u - u') + (Ai - A) / L) + (1 - 6ab) s I L, (2.8) 

and we have three variational parameters A, Ai, and s. The variables A, Ai represent two 
values of \{uj ^ 0) and Xiu = 0) instead of the general function A(co') introduced in Ref. 
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^ (which involves an infinity of variational parameters). The variable s represents a "spin 
jlass" type variable which loosely speaking is a measure of "freezing" . 
The variational Hamiltonian now becomes 



hn — Jo du J2a=l 



+ jr EabPab Jo du Jq du' Kaiu) " Rfe(M'), 



(2.9) 



with 



/ 



Pab 



\ 



V 



(2.10) 



Ai — A —s ■ ■ ■ —s 
—s Ai — A ' ■ . : 
: ' • . ' • . —s 

s ■ ■ ■ — s Ai — A 

which reduces to the EM variational Hamiltonian if pab = (i.e. if Ai = A and s = ). 
For now we consider a replica symmetric parametrization. We will discuss a possible replica 
symmetry breaking parametrization later on. Using this parametrization of hn our task is 
to calculate the free energy from equation ( |2.6|) . This is achieved by first writing down the 
propagator associated with phn. 



Ma;2 + ^)l-q(^))}~^' 



(2.11) 



with 



qabi^^) = / du qab{u) exp{-iuju) . 



(2.12) 



For the function qab{u) given by Eq. (|2.8|) we find: 



qab{^^) = -Sab ((A - /U) + (Ai - A) S^fl) + (1 - Sab) S S^^Q, 



(2.13) 



and thus 



which gives after inverting an n x n matrix and taking the limit n ^ 0, 



(2.14) 
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PGabiu; = 0) = .^'^\ Sab + ' ,, (1 - 6ab), (2.15) 
(Ai + s)^ (Ai + s)^ 

PGab{u;j^O) = ^J^^ 6ab. (2.16) 

Since the interval on which u is defined is finite [0 < u < L ), the "frequencies" u are 
discrete and satisfy 

a;^ = ^m, m = 0,±l,±2,.... (2.17) 

We can now use the fact that 

{Ra{u) ■ Rbiu')) = d g^biu -u') = yY. e-^^^^-^'^G^fe (o;) , (2.18) 

to obtain an expression for the correlation function of interest and for the free energy. For 
the correlation function we obtain 



Cii) 



1 1 " 

= - (((R(£)-R(0))^))^ = - 5: ((R„(£)-R,(0))^ 



For the free energy we find from Eq. ( |2.6| ) 

n (F) 1 " 1 

— — = const. + -(/i - A) ^ ^GaaV) - - ^PabGab{i^ = 0) 

a=l tt) ab 

-4Tytr In G(uj) 
2/3 V 

/"^^ f dk ^ 



2d 7o J (2vr) 
/ k2 



exp 



2L 



[Caaiuj) + Gbbiuj) - 2e~''''Gabico)] . (2.20) 



We now use the formula (see e.g. Gradshteyn and Ryzhik |jT^, Eq. 1.445.2) 

ly- _ 1 cosh(a(l - 2^/L)) ^ _ L [Y ..... 

Mu;^ + X~ 2VMX sinh(a) ' ""2Vm' " " ^ ^ 

to calculate the correlation function and the free energy in the limit n — 0: For the corre- 
lation function we obtain 
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Cm = ^— Icoth^ cosh(L^(l-2£/L)/2) ^ 

/^TmA I^""" 2 V M sinh(LyA7M /2) ' ' 



and for the free energy 



^ const. + coth #,/A + ^ ^ ^1+2^ - i 



+ |lnsinh#./5+^lnfl + f) + ^ln^ ^ 



L "^'^"^^^ 2 V A/ ' 2L V ^1/ 2L A 2L Ai+s 

t9 [L J r _dk 

2d JO J (27r) 



-^lodz Jj^ [exp(-k2ai) - exp(-k2a2)] , (2.23) 



with 

1 I [\ cosh{L^X/M{l-2z/L) /2)\ 

ai = ^= coth-W ^- ; , (2.24) 

2I3VM\\ 2VM sinML^A/M /2) / 

1 / 1 1\ 1 L [T , , 

a^= [ + ^coth-W— . 2.25 

Some of the details of the calculation are given in the Appendix. The constant term does 
not depend on the variational parameters. So far the calculation has been exact but now 
we are interested in the large L limit. Before we proceed it will be instructive to pause to 
review the calculation of EM who have chosen p = 0, i.e. s = and Ai = A (recall that in 
their notation A oc g^). They also take /i = 0. In that case the free energy simplifies to give: 

^^"^^ - const. - -4^coth-^./^+ Tlnsinh# 



^ 4VMA 2 V L "-'^"--^ 2 V A/ 

- ^ lo^' dz J [exp(-k2ai) - exp(-k2a2)] , (2.26) 

with Oi still given by Eq. (|2.24|) and 

02 = ^= coth -\ —. (2.27) 

2PVMX 2\ M ^ ^ 

We also noticed that since ai is symmetric about the point z = L/2 we have limited the 

z-integration up to L/2 and multiplied the integral by 2. We can now take the limit of large 

L . It is at the point z=L/2 that the integrand vanishes for large L. In this limit we find 

(upon dropping the constant): 

(3{F) 1 f3^g , f dk 



Ld 4 V M d Jo J {27tY 
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poo r dr^ r 1 

J dz J — — [exp(-k2ai) - exp(-k2a2)J , (2.28) 



with 



ai 



1 — exp 



^\/Vm))' 



1 



02 



(2.29) 



Notice that the factor of 2 in front of the integral due to the aforementioned symmetry was 
missed in Ref. This is of no importance since it just renormahzes the strength of the 
disorder. The integral over k can now be done to yield: 

d/2 



I3{F) _l X f^^g(p^fM\\ 



Ld 



d \ 2n J 



dz 



1 



d/2 



- 1 



(2.30) 



^1-exp (-z^A/m)^ 

At this point we realize that the 2;-integral is infrared divergent for any dimension d > 2. We 
can trace this back to the short distance singularity of the Dirac delta function correlation. 
We thus replace the delta function by a regularized form: 

1 f 



exp 



(2.31) 



where ^ is small (we can think of it as the intrinsic diameter of the polymer thread) . Using the 
representation given by Eq. ( ^.41 ) for the right hand side and carrying out the y-integration 
yields: 



Ld 



- — 

d Jo 



dz 



exp(— k^ai) — exp(— k^a2) 



(2.32) 



and the integrals are now properly regularized to yield a finite expression. To find the 
optimum variational parameter A, we take the derivative of the above expression with respect 
to y/X : 

, 2/5(7 IM 



dX 



X 



(1 — exp 



-r — r 



exp (— t)) exp(— k^ai) — exp(— k^a2) 



(2.33) 



and the variable has been rescaled hy z ^ ryM/X. The k-integration can now be done 
to yield: 
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1- l^^2 {{1- + Af'^' ~ (1 + A)^/2+i j ' ^^-^^^ 



with A = At this point we see that the integral is finite for d < 4 even in the 

hmit A — s> (which follows from ^ 0). Let us denote the integral in this limit by 1^- 







so, Eq. ( p.34|) becomes 



d-4 



Id A— (2.36) 



(27r)'^/2 

Unfortunately EM did not realize that this integral is negative for d = 3 (and also d = 2). 
The indefinite integral can be carried out analytically (e.g. using Mathematica) in d = 3 to 
give 



Is = lim^^oo Ht) - lim^^o /(^) = -§(!- In 2) -0.20457. (2.37) 

(In d = 2, one has I2{t) = t/ (e^ — 1), and I2 = —1). From here it follows that Eq.( p.36| ) has 
no solution for A ! This is a very important observation. Notice that all fractional powers 
of A are always to be taken as positive. For example in the integral 

= 23437i(2;3VMA)'"2+i, (2.38) 

which is part of the result derived above, in the integrand is positive, and so must be 
the result of the integration since the integrand is positive definite. There is no way to argue 
that A^/^ can be taken as the negative square root of V^. To elucidate further the fact that 
there is no value of A which extremize the variational free energy we return to Eq. ( p. 32 ) of 
the free energy and carry out the k-integration to find for (i = 3 

3L ~ iVM 3(27r)3/2 

I ( (1 + A - exp(-r))3/2 - (1 + A)3/2 ) ' (2-^^) 
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with A = i'^dpVM\/2. A gain the integral can be done analytically (Mathematica) and we 
find 

substituting for A one obtains 

We see that the divergent term (as ^ ^ 0) is independent of A. We also see that the free 
energy is a monotonically increasing function of A and thus has no extrema as a function 
of it. Derivative of the last expression with respect to agrees with our previous result. 
Thus we see that the one parameter variational Hamiltonian does not yield a meaningful 
result. 

Let us now return to the more general expression for the variational free energy given in 



Eq. ( p.23|) . Before we consider the large L limit we can draw a general conclusion. Let us 
calculate the derivative of the fee energy with respect to Ai and s: 



\3 



2L(Ai + s 
' da2 \ , f dk 



dz / -^k2exp(-k2a2), (2.42) 



\d\i) 2d Jo J {27tY 
s (2/x — Ai — s) f da2\ jS'^g r dk 



2L (Ai + sY \ds J 2d Jo J {2n 



Since 



we find that upon adding the two equations we get 

2L (Ai + s) 

which implies 

Ai + s = fi. (2.46) 
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This is an important general result. Substituting this result in Eq. (|2.42| ) we find 

f^9 r f dk 2 ^ 27T(3gL 



s = / - — -;k^exp(-k^a2) = 7 , , (2.47) 



with 



02 = — --- + ^coth-W— , 2.48 

which gives a relation between s and A. However as we will see in a moment, only the 
combination Ai + s = yU, enters in the equation for A. 

Returning to Eq. (p.23|) we see that upon taking the derivative with respect to -\/A the 



only dependence on s and Ai is through the combination Ai + s in a2. It is simpler to take 
the limit of large L and to write up the resulting equation for A up to exponentially small 
terms in L : 



X (exp (— k^ai) (1 — exp (— r) — rexp (— t)) 

-exp(-k2a2)(l-4./M7A/L)), (2.49) 



with 

1 - e"^ 1 11 , , 

«l = «2 = ^= + 777 77TT- 2.50 

If /i is finite, one can proceed with expanding exp (— k^a2) in powers of 1/L as will be done 
later. However if one attempts to take the limit /i — we see immediately a potential 
problem because of the term l/(/5L/i) in 02. If we carry out the k-integration we find 

2(7/3'^/2+2M'^/4+i ^ , 
= (2^^ 



(2.51) 



where we have omitted subleading terms in 1/L. As /i — » for fixed large L, the last 
term in the integral vanishes (as A — /i remains finite for g ^ 0). The integral over r no 
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longer converges for large L, but is rather proportional to L. To leading order we get (by 
subtracting and adding 1 to the integrand) 



d+2 



which gives 



4 
2-d 



(2^^ ) ■ ^^-^^^ 

We see that the borderline dimension appears tohe d = 2. Indeed from Eq. ( p.22|) it follows 
that for large L 

m ^ (l - eM-^\fm)) , (2.54) 

and thus the radius of gyration satisfies 

^ A"^/^ ~ (^L)-^, (2.55) 



which agrees perfectly with equation (|1.6|) for d < 2. To see what happens for d > 2 we can 
easily show that in the limit fi the free energy becomes of the form 

^ rf . 1 



P{F) ^ 1 /A _ ^ ( /3VMX \ 
Ld A\M d [ 27r ) V A 7o (1 + A - exp (-r))'/' 



dr--- -— ^, (2.56) 



where again we reg ularized with A = i'^dp^fM\/2. This gives 

m.,on=t. + i,/A_£!^fi!^y'\, (2.57) 

Ld 4 V M 2ci V 27r y ' ^ ^ 

and when using A ~ d'^ (3~'^ we obtain 

,5(F) = co„st.xL + — (2.58) 



which coincides with Eq. ( |1.4| ) and shows that for(i>2,F^— ooasi?^0 and there is 
always collapse. Thus we see that in the limit of /i we recover the annealed result from 
the replica calculation as expected. 

If on the other hand /i is finite, we can expand exp (— k^a2) in powers of 1/L and we find 
to leading order in L: 



The last two terms are also 0(1) although they originated from seemingly 1/L terms, since 
we obtain a factor of L from the range of integration over r. Evaluating the integrals we 
find 

For small g we can solve this equation in powers of g. Defining a dimensionless constant 

^ = ^;^p7^ ' (2-61) 



we cast the Eq. ( p.60|) in the form 



h = l + g2{h + ^ + '^-^h)h"\ (2.62) 

with h = X/fi. To second order in g we find: 

X/fi=l + 2{h + 2)g + (Jrf + 2) (d(/rf + 2) + 2(ci + 2)) ^2 ^ . . . . (2.63) 

Thus as g increases from 0, A is an increasing function of g starting from an initial value 
of /X. However a numerical solution of equation ( p.60|) (for d = ?>) reveals that the solution 
becomes ill behaved as A becomes of magnitude ~ 2/i. This happens for ^ ~ 1/ [2^/ ^{1^ + 2)). 
The reason for this is as will become evident in the next section is that the replica symmetric 
solution becomes invalid at this point and has to be replaced by a replica symmetry breaking 
solution. This will become clear in the next section where we will find the correct solution 



for larger values of g. It is also clear from Eq. ( p.61| ) that for fixed g as fi 0, g becomes 
large and we will be in the region when RSB is to be used. Thus the range of applicability 
of the replica symmetric solution is minimal for a small value of fi. 
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The rest of the section can be skipped on first reading of the paper and the interested 
reader might continue directly to the next section discussing the RSB solution. For com- 
pleteness we display here the form Eq. (|2.59| ) takes for a general correlation of the disorder 



defined in Eq.( |1.14D . We can use the representation given in Eq. ( p.4[ ) to obtain 



oo 



A — /i = —A/3g\l —I dr 



A Jo 

x(f'( (1 - exp (-r) - r exp (-r)) - / ' 



where we defined 

/(a) = / dy /(y Vrf) / exp(-zk ■ y) exp {-^) 

= TmIo^dxx'^/'-'e-^f (^), (2.65) 

and the primes stand for derivatives of /, which can be obtained from the first line of Eq. 
( p.65| ) by taking the derivative with respect to a under the integral sign. 

At this point we would like to discuss the more complete variational scheme that we used 
in Refs. |^,|^ and show that all our conclusions concerning the limit /i ^ follows from that 
scheme as well. Since the notation there was different we will translate the equations to the 
present notation but we will not rederive them here. What we did there was to consider a 
variational scheme in which we allowed the variable A to depend on u and we extremized the 
free energy with respect to each variable X{uj). The propagator G{uj) defined in Eq.( [^.14[ ) 
now becomes 

Gabiuj) = {{Muj' + X{u) + (Ai - A(0) + s) 5^,o)l - s 5^,o)}^^' • (2.66) 
We have found that the relation Xi + s = fi still holds and s and A(cij) satisfy the equations 

s = -2(3gLf'{2a2), (2.67) 

X{uj)-fi = -s-2/3g dz {l-e''^')f'{2ai{z)), u^O, (2.68) 

Jo 
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with 

1 ^ 1 - e-'^'' 



For a regularized delta function correlation we have 

^ '^""^ " ~2(27r)'^/2 (de/2 + a)'^/2+i- (2-71) 
In the limit /i — 0, we observe that s — 0, and there is no longer a cancelation of the 
contributions linear in L between the two terms on the right hand side of Eq.( |2.68| ). Instead 
we get 

AM^-2^.L/-(A^^^^^^). (2.72) 

which yields an u- independent solution that for the delta correlation becomes 

This result exactly coincides with Eq.( [^.52[ ) derived previously. 

III. REPLICA SYMMETRY BREAKING 

In the previous section we have seen that the replica symmetric solution becomes invalid 
for fixed amount of disorder and small harmonic constant fi. In this section we show the 
emergence of a different solution of the variational equation which is more adequate for our 
problem. But in order to take advantage of such a solution we must use a more general 
variational scheme. Returning to Eqs. (|2.8|) - (|2rT0|) , we have extended the parametrization 
of the matrix pab in ( |2.10| ) to allow for one-step RSB by having two off-diagonal parameters 
So [x < Xc) and 51(2; > Xc) together with a breaking point Xc (0 < < 1) . Here x is 
Parisi's replica index. For details of Parisi's RSB scheme see reviews of spin glass theory 
P^3|^T7[. Thus our variational scheme includes now 5 parameters. A one step breaking is 
sufficient for the case of short range correlations of the random potential @,|l^ . 
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We were able to calculate analytically the free energy with the new parameters. Here we 
display the final result, the details given in the Appendix : 



P (F) (/i - A) 



A 



Ld 



So 



coth —\ — + —r 



2 V M 2L \Xc{Xi + si - S) 



1 



A 



1 



(Ai + si-S)2 A 
i ( 1 - - ) In ( 1 - 



- + -Insinh-W — + — In 1 + 



2L 
2d 



L 



2 V M 2L 



Sl 



Ai 



Ai + Sl 



1 



So 



1 l^Ai 
2L A 2L Ai + Sl - S 



dz 



dk 

{2Tiy 



exp(— k^ai) — Xcexp(— k^a2i) 



Ai + Sl 



Xc) exp(-k^a26) 



(3.1) 



We introduced the notation 



E = Xc{si^So) 



(3.2) 



the variable ai is still given by Eq. (|2.24| ) and we defined 
1/1 1 / 1 \ 1 1 



a2i 



j3L Vxc Ai + Sl — S 
1/1 r 



1 M ^ 
~^^/A7+si \J ' 2/3yMA 2 V M 

1 , L fx 

coth ■ 



1 1 ^ / A , , 

coth— W— , (3.3) 



(3.4) 



/?LVAi + si A/ 2/5VMA 2\ M 

From the free energy we are able to get the following five relations (everywhere we eliminated 
Sl in favor of S) 



Ai + So - (1 - l/xc)S = /X, 



(3.5) 



which replaces the relation Ai + s = /i established above for the replica symmetric solution. 

(3Lg f dh 



So 



d J {2tt) 



(3.6) 



^ — — "^LXf [ ~~ — TT-I 

d J {27rY 



exp(-k^a26) - exp(-k^a2i) 



(3.7) 



^ ln(l I ^]-^'Lx ^ f 
/X + S [ fij d 'fi{fi + J:)J i2nY 

H — -{Lxc) 



d 



dk 



{27iy 



k^exp(-k^a2«) 



exp(-k^a2i) - exp(-k^a26) 



(3.8) 
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^ d ^ \ Jo J {2nY 

X (^exp (^—li^aij (1 — exp (— r) — rexp (— r)) 



(-Xeexp (-k^as/) - (1 - Xe) exp (-k^asf,)) (^1 - A,/M/X/Ljj , (3.9) 
where we defined 

1 - e"^ 

ai = = (3.10) 

_ 1 1 1 

" 2/5v/MA ~ ^ + /3L(/. + E) ' ^^-^^^ 
1 111 1 /I 1 \ , , 



2(3 VmX (3LX (3fiLxc /5(/i + S)VL Lx^ 
We have simphfied some expressions assuming large L and dropped a term of order 1/L in 
Eq. (^). 

If we denote by i/c = Lxc we reahze that equations ( |3.7|) and ( p.8|) can be solved for S 
and Uc of 0(1) with respect to L. These equations are similar for those of a classical particle 
in a random potential, |]TH[ except for the variable A which does not appear there. (One can 
recover the equations for the classical particle by taking the limit M — oo with L fixed. 
One needs to replace (3L with j3 for a particle. This limit is not meaningful for a polymer.) 
For small fi we can have an approximate analytical solution: 

-(/^VmA) V^W^, (3.13) 



{27tY/' 



Vc ^ Lx, = ^W-(27r)'^/^(/?%/MA)-'^/^^^h^, (3.14) 



So 



const. X ^7(2-'^)/4/3L(/5yMA)-"('^+2)/V'/'+' |ln /^|^'+'^/' • (3.15) 



An analysis of the equations (expanding in power series in S) shows that this solution is 
valid as long as the condition 

2/5^(1^^1)"'*'%"""' SI (3.16) 
is satisfied. This inequality can also be expressed in the form 

g(d + 2)i-j >1, (3.17) 
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where g has been defined in Eq. ( |2.61|) . When the equahty holds we have S = and Xc = 
(4 + d)v/MA/(2/iL). This can also be verified by using this condition at the equality point 



in the above solutions for S and Xc and we see that indeed S ~ 0{fi), and Xc ~ V MX/ {fiL). 
Solving the equality condition given by Eq. ( |3.17| ) together with Eq. ( p.62| ) gives 



h = l + y ~ 1-85 for d = 3, (3.18) 
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{d + 2)/i(^+^)/^) ^ ~ 0.068 for d = 3, (3.19) 



in agreement with our numerical solution of Eq. (|3.17D which broke down for h ^ 2 ior d = 3 



So the point when the replica symmetric solution has to be replaced by the RSB solution is 
just below the point that the RS solution becomes ill behaved. 

If on the other hand /i is small but fixed we can use the solution we have obtained above 
in the equation for A, in the limit of large L. We obtain 

2(3g [m ^ f dk 
A — /i = — — y— / dr 



d V A Jo J i2ny 
(exp k^ao(l — c^^)) (1 — exp (— r) — rexp (— r)) — exp k^ao 

(3gyc f dk ^ ( ( .2 ( , ^ \\ /,2 
k exp — k ao + -. =;r — exp ( — k ao 



d J {27tY V V V /?Z/c/i(/i + S), 

A(3g [m r dk . / , 2 \ 



(27r)' 

d + S ' a) 7 {2tiY' 



+ -A^-\]I T^^k^ exp (-k^ao) , (3.20) 



where we defined 



ao = l/(2/?yMA). (3.21) 

We can check that for S = it reduces to the replica symmetric equation. Carrying out the 
integrals we get 

2(7/5'^/2+2M'^/4+^ 

A — fJL = - 



(27r)^/2 

'k [X 

2\My y ' y^fi{fi + S) y ) ' 2 fi + E 

24 



xU^ + 2 + ^Jt7 1- 1 + — 77-^ ]+—. Ai, (3.22) 



and we have to substitute for S and Uc (which are functions of A) from Eqs. ( |3.13|) and 
( p.l4| ) respectively. If we now consider the case of strong disorder we can neglect n relative 
to A and S and the above equation simplifies to give 

Substituting for S and Uc we find 

_ 2g(3'^l'^'M'^l'^' ( 2-d 
(27r)'^/2 2 



^ Vf/^ (2^)i'^/5-(4+<^)/4M-(^+'^)/«A(^-'^)/«V^ ( 1 + 4^] I • (3-24) 
Let us seek a solution of the form 

A = C^/^^-'^) {(3'm) ^'^'^^^'-"^ gm-d)_ (3.25) 

Substituting in Eq. (^.241) we obtain a quadratic equation for C and to leading order as 



^ we find: 

^ = {2Tifd/iA-d) ^) (9 |ln/x|) . (3.26) 

Using this result inside the parenthesis in Eq. (|3.24|) we see that we get + 2 + |ln/i| 



. This shows that we were justified a posteriori in neglecting the constant terms. It also 
shows that the negative constant Id of EM (see Eq. |2.36|) has been replaced by the term 



|ln fi\ . From this final result we obtain the radius of gyration 

-l/(4-d) 



^ {^P^g lnVy'^''-'\ (3.27) 

This is the main result of the paper. It recovers the EM result but with their constant 
Id being replaced by 21nV^ as has been argued by Gates and Ball [@]. Note that we have 
replaced M in favor of the bond step b. 

Substituting the result ( p.26|) obtained for A in Eqs. (|3.13| and (|3.14| ) we find 
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S = {-^.gP^'^-'^'^M^^^^y^ |ln/i|) (3.28) 
Vc = Lx^ = {^g'P'^'M'^ \\n^^tY^^'~'^ • (3-29) 

The second equation is important since Xc can not exceed 1 (Parisi's variable x must satisfy 
< X < 1 [|1^]) . For 2 < (i < 4 we see that Xc actually decreases when /x becomes small 
so there is no problem. Also for d = 2 there is no problem since for large enough L, Xc 
is also within range. On the other hand when d < 2, Xc increases when /i becomes small 
(or equivalently V becomes large) and eventually will exceed 1. For example for d = 1 we 
see that this happens for |lnyu| ~ g^L^, which corresponds to an extremely large volume V 
~ exp((7^L^) when L is large. For V > V we revert to the annealed result, which for d < 2 
predict R ~ (LgY^^'^^'^^ as was shown in the last section. In the large L limit this again 
leads to a fully collapsed polymer. 

We have also verified that to leading order the free energy is given by Eq. (|1.1CI| ) (there 
is a subleading term of the form Lg/R^'"^ that can be neglected). It is interesting that the 
condition Xc < 1 that we have applied above has a physical significance [0]. The attractive 
term in the free energy is (see Eq. ( |1.10| )) of the form —L^g \n.V/ R^. This represents (up 
to a sign) the binding energy of the chain. In order that the polymer will be confined to 
a small single region of size R as given above in Eq. ( p. 271 ), the binding energy should not 
exceed the translational entropy ~ \n.V . The condition 

\nV < L^g\nV/R'^ (3.30) 

is equivalent (up to some irrelevant constants) to the condition Xc < 1 as can be verified by 
using the result in Eq. ( |330|) . 

IV. CONCLUSIONS 

We have considered the problem of a polymer ( a Gaussian chain) in a quenched disor- 
dered medium. The problem maps also to a quantum particle in a random potential, and in 
the presence of an additional confining harmonic force (of spring constant fi ) it maps also to 
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the problem of a flux line in a cage potential and random columnar disorder. We carried out 
a replica calculation in the presence of a confining harmonic force, and succeeded to "im- 
prove" the previous results of EM |jl|] , in the sense that the (unphysical) constant is replaced 
by In V in the equation for the variational parameter A and hence also in the dependence 
of the radius of gyration on the strength of the disorder. Of course our calculation does 
not diminish the accomplishments of EM who pioneered the use of the variational method 
in the context of the replica calculation and for the first time obtained the correct scaling 
exponent for the dependence of the radius of gyration on the disorder for finite systems. In 
the infinite volume limit the chain collapses since it can find a very deep potential minimum 
somewhere which can accommodate it. For 2 < d < 4 the chain is localized in the sense 
that even in the large V limit two long chains introduced into the system will find the 
same small neighborhood to occupy (with a probability approaching 1 for large L). This is 
a consequence of the off diagonal spin-glass order parameter we introduced that measures 
overlap between different replicas. It is comforting to find out that the replica calculation 
can reproduce all the physical arguments introduced so cleverly by CB [Q. 
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APPENDIX: 



Here we give some of the intermediate steps leading to Eqs. (|2.6|) and (|2.2CI|) . To evaluate 
the expectation value of the last term in Hn it is most useful to write: 

(expik- (R,(u)-R,(M')))ft„ = / VKi-'-VKnexp [ V,(i;)-R,(i;) 
-\ E dv dv' R^iv) g-J{v - v') R,{v')\ 
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(^1 PRi • ■ ■ PR„ exp (^-i £ dv dv' R,(t;) g^,\v - v') R,(i;') 



where 



exp ( --k^(5(„„(0) + gbb{0) - 25(„fe(u - u] 



V,iv) = I k {b^J>{v -u)- - n')). 



Next we show how to evaluate other contribution to the free energy: 

L h [\ \ Ai + 2s 

coth — W — ~ V + 



Also 



A(Ai + 2s) 1 



- ^EPabPGabiio = 0) = -^Tr p G(0) 



ab 



2n 



2(Ai + s)2 2' 



and 



- Vtr ln/3 GM = ^Etr In (r^G-i(^; 



E ln(Mo;2 + A) - - In A + -tr In (/^-^Q-^ 



but 



Eln(M^2 + A) = ^ In |^2sinh 



n const. 



(Al) 



(A2) 



(A3) 



(A4) 



(A5) 



(A6) 



(see e.g. |]T^ p. 44, Eq. 1.431.2). The constant term (which is infinite) is eliminated by the 
normalization of the functional integral, and in any case does not depend on A. Also 



tr In [(3-^G-\0)) = tr In 



Al — s 
-s Al 



V 



n In Al + n In 1 H — n 

* Xj Xi + s 



■ — s Al 
+ o{n^). 



(A7) 
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For the case of 1-step RSB we have to calculate the propagator by inverting Parisi type 
matrices. It is helpful to use formulas found in an Appendix of [0. We find 



a;c(Ai + Si - S) V a;c/Ai + si (Ai + si - S)^ 
■So 

(Ai + si-S)^ 



= 0, x) = , ^° , X < Xe, (A9) 



PG{u; = 0,x)= ^ — - - + (^ ^ v^2 ' ^ > (^^0) 

Xc(Ai + Si-S) XcAi + Si (Ai + si - S)^ 

^ 0) = ^J.^^ Sab , (All) 

where x is Parisi's index on the interval [0,1]. 

Next we show how various other contribution to the free energy become in the RSB case: 

- ^T.Pabf3Gab{uJ = 0) = -i^Tr p G(0) 

ab 

= ^PGaaiuJ = 0) - i (A13) 



and finally 



2 ' ' 2' 



■^Etr ln/5 GH = i-^tr In (r^G-^(..; 
:^5:in(Mo.2 + A)-^lnA + ^tr In (^^0-^(0) 

lnl2sinh -\[^]+-\n^ + -ln(l + '''^ 



2VM/2A2V Ai 

The coefficient ai in the exponential is given as before by 

ai = jY.Gaa{oo){l-e'^^), (A15) 

and a2 becomes 
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a2{x) = \ {Gaai.^ = 0) - G{U} = 0, X)) + I J2w^0 G'aa(^) 



a2l, X < Xc 



026, X > X,. 



(A16) 
(A17) 



and the explicit expressions for oi, a2i and 02?, are given in Eqs. ( p.24| ), (|3T 
respectively. Notice also that 



and 



exp(— k^a2(a:)) = — / rfx exp(— k^a2(x)) 

= -Xcexp{-k'^a2i) - (1 - x^) exp(-k^a26) 



(A18) 
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